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摘 要 : 本 文通 过 使 用 各 向 同性 消 流 燃烧 的 可 压缩 DNS By 


发 现 忽略 密度 脉动 和 三 阶 脉动 关联 量 是 可 行 的 ， 可 以 有 效 地 简 


重点 实验 室 ， 浙 江 杭 州 310027) 
居 库 对 二 阶 矩 (SOMD) 燃 烧 模 型 进行 了 先 验 性 检验 和 评价 。 通 过 检验 ， 
化 反应 速率 方程 。 通 过 对 关联 量 输 运 方程 中 的 精确 项 的 检验 ， 发 


现 化 学 反应 项 的 贡献 最 大 ， 需 进行 准确 的 封闭 。 通 过 对 模 化 的 输 运 方程 的 检验 ， 发 现 部 分 精确 项 的 模 化 效果 不 理想 ， 但 总 体 来 


Ui, SOM 燃烧 模型 的 思想 和 方法 是 可 行 的 ， 有 待 进一步 的 发 


关键 词 : 二 阶 窍 ， 消 流 燃 烧 模型 ， 直 接 数值 模拟 ; 先 验 性 验证 
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A-priori Validation of a Second-order Moment Combustion model 
BAIYun LUO Kun YANG Jian-Shan QIU Kun-Zan FAN Jian-Ren 
(State Key Laboratory of Clean Energy Utilization, Zhejiang University, Hangzhou, 310027, China) 


Abstract: A turbulent combustion model, Second-order Moment (SOM) combustion model is validated a-priori 


using the DNS database of combustion in isotropic turbulence. With the validation it is found that neglecting the 


density fluctuation and third-order fluctuation correlation terms can simplify the reaction rate equation efficiently. By 


comparing the exact terms in the correlation terms transport equations, the reaction terms play a dominant role which 


need to be modeled more accurately. Though some predicted terms are not in good agreement with the DNS values, 


the SOM closure method can be a feasible method for turbulent combustion modeling and it is necessary to be 


developed. 
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淇 流 燃 烧 广泛 应 用 在 能 源 、 冶 金 、 化 工 、 交 通 
等 工业 领域 ， 满 流 燃烧 的 研究 ， 有 助 于 提升 对 燃烧 
的 认识 和 应 用 。 随 着 汕 流 研究 和 计算 机 技术 的 发 
RANS, LES 和 DNS 成 为 了 研究 满 流 燃烧 的 重要 途 
径 和 方法 ， 并 且 各 种 燃烧 模型 被 提出 和 应 用 。 
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项 的 封闭 ,学 者 们 提出 了 许多 不 同 的 满 流 燃烧 模型 。 
如 EBU-Arrhenius 模型 中 概率 密度 函数 (PDF) 模 
型 外 ， 小 火焰 模型 站 ， 输 运 PDF BAO, ZEKE 
(CMC) 模 型 "I 等。 但是， 上述 的 燃烧 模型 ， 都 是 针 
对 不 同 的 燃烧 模式 提出 的 ， 难 以 同时 预测 预 混 火 焰 
和 扩散 火焰 。 同 时 ，DNS 计算 表明 09， 预 混 火 焰 和 


型 


扩散 火焰 会 同时 存在 ， 这 给 上 述 燃烧 模型 提出 了 挑 
战 ， 能 够 同时 预测 不 同 燃 烧 模式 的 模型 成 为 了 燃烧 
模型 发 展 的 一 个 重要 难题 。 近 年 来 ， 由 周 力 行 哺 ] 提 
出 的 二 阶 和 矩 燃烧 模型 (Second-order Moment, SOM) 
将 尝试 解决 这 一 问题 .SOM 燃烧 模型 是 基于 反应 速 
率 的 Arrhenius 方程 ,将 温度 的 非 线 性 指数 函数 用 一 
个 反应 率 系数 开 表示 , 然后 由 基本 输 运 方程 推导 得 
到 脉动 关联 量 的 传输 方程 ， 进 而 采用 类 似 满 流 封闭 
的 方法 对 其 输 运 方程 进行 封闭 。 理论 上 , SOM 燃烧 
模型 对 燃烧 模式 没有 限制 , 因此 , 可 以 对 预 混 火 焰 、 
扩散 火焰 ， 以 及 部 分 预 混 火焰 进行 很 好 地 预测 ， 同 
时 ， 对 于 着 火 和 熄火 等 强烈 化 学 反应 现象 也 可 以 进 
行 捕捉 。 周 力行 等 ("的 研究 中 , 认为 EBU-Arrhenius 


模型 以 及 简化 PDF 模型 在 理论 上 存在 缺陷 , 模拟 结 
果 与 实验 值 出 入 较 大 。 而 PDF 输 运 模型 、 小 火焰 模 
型 以 及 CMC 模型 ， 模 型 较为 复杂 ， 计 算 量 较 大 ， 
而 且 ， 不 能 同时 预测 扩散 和 预 混 火 焰 。 因 此 ， 相 对 
简单 的 、 基 于 反应 速率 封闭 的 SOM 燃烧 模型 成 为 
了 有 潜力 的 消 流 燃烧 模型 ， 有 望 在 工程 中 模拟 消 流 
燃烧 ， 并 在 相关 文献 中 得 到 应 用 。 

然而 ,， SOM 燃烧 模型 存在 一 些 假设 , 包括 忽略 
密度 脉动 和 忽略 三 阶 脉 动 关联 量 ， 并 且 模 化 结果 需 
要 检验 。 因 此 ， 本 文 使 用 可 压缩 DNS 数据 库 对 模 
型 进行 了 先 验 性 检验 ， 评 估 了 SOM 燃烧 模型 。 


1 SOM 燃烧 模型 


对 于 双 组 份 的 单 步 一 阶 总 包 反 应 ， 其 瞬时 反应 
速率 的 Arrhenius 形式 为 : 


w, = 了 DT 了 了 exp(—E/RT) (1) 


其 中 ， 下 标 应 表示 燃料 ，ox 表示 氧化 剂 ，Yi 
和 Yo 分 别 表示 两 种 反应 物 的 质量 分 数 ,B 为 方程 系 
数 ，E 为 活化 能 ，R 为 通用 气体 常数 。 

对 方程 (1) 进 行 Reynolds 平均 并 : 
和 三 阶 脉动 关联 量 得 到 时 均 反 应 率 为 : 


忽略 密度 脉动 


w, = Bp (Y, Y K+ KY Y, + Y, KY, + Y.KY,) (2) 
程 (2) 中 的 浓度 -浓度 关联 量 YY 和 浓度 -反应 率 系 


其 中 ，K = exp(-E/R7)， 
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数 关联 量 KY' 、 


关联 量 的 精确 输 运 方程 通过 组 分 方程 和 能 量 
方程 推导 得 出 ， 其 中 浓度 -浓度 关联 量 精确 方程 为 : 


KY! 通过 求解 其 输 运 方程 得 到 。 


OOY, Ya PUY Y Ou Yp Ya 
ôt Ox Ox. 
# ei 
Ol, oy. a RAA 
-| pu,’ "+ pu Y, e |+ pb 5 (3) 
Ox, ox, Ze, 


式 中 ， 右 端 项 的 第 1 AMT LL, 3 2 项 
为 平均 浓度 梯度 引起 的 产生 项 ， 第 3 项 是 层 流 扩散 
项 ， 第 4 项 为 分 子 耗 散 项 ， 第 5 项 为 与 化 学 反应 有 
关 的 源 项 ， 这 些 项 都 需要 进行 封闭 。 其 封闭 后 为 : 


模 化 方程 右 端 项 第 
生 项 ， 第 3 项 为 耗 散 项 ， 


’ Cer Ce 和 Co 是 经 验 系数 ， Tr FAT, 分 


流 扩 散 时 间 和 化 学 反应 时 间 ， 定 义 为 : 


4 
RT 


同 理 得 到 浓度 -反应 率 系 数 关联 量 的 精确 方程 : 


1 项 为 扩散 项 ， 第 2 项 为 产 
第 4 项 为 化 学 反应 项 ， 式 


PAW Ania 


apK'Y' pu KY) pu KY' D ô En 
ôt Ox, ox, ox, Ox, 
ay EPET aK ðY’ 6K! 
=| pu K'— + pa 一 |-2pD 一 一 (6) 
ax, ox, ox, Ox 
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与 浓度 -浓度 方程 一 样 , 右 端 项 的 第 1 项 和 2 项 
为 扩散 项 ， 第 3 项 为 产生 项 ， 第 4 项 为 耗 散 项 ， 第 
5 至 第 8 项 为 化 学 反应 源 项 。 采 用 类 似 的 方法 封闭 ， 
模 化 方程 为 : 


nA 
ae ore (ae a [4 ra | 


类 似 的， 方程 右 端 项 第 1 项 为 扩散 项 ， 第 2 项 
为 产生 项 , 第 3 项 为 耗 散 项 , 第 4 项 为 化 学 反应 项 。 

上 述 就 是 二 阶 矩 燃烧 模 MAEA 思想 ， 需 要 对 
其 假设 和 模 化 方法 进行 检验 。 王 方 等 "4 使 用 村 道内 
不 可 压 等 温 反应 流动 的 直接 模拟 对 SOM 燃烧 模型 


F 


进行 了 检验 ， 本 文通 过 浙江 大 学 能 源 清 洁 利 用 国家 
点 实验 室 提供 的 可 压缩 DNS 数据 库 对 该 思想 方 
法 进行 先 验 性 检验 ， 并 作出 评估 。 


2 DNS 数据 库 的 描述 


本 文 使 用 的 DNS 数据 库 是 王 海 欧 等 [5 采用 单 
步 不 可 逆反 应 机 理 ， 直 接 数值 模拟 计算 了 各 向 同性 
MAPAS R E. 计算 区 域 沿 x, > 和 = 方向 都 是 
L=2x mm， 网 格 数 为 128; ， 网 格 大 小 为 494m ， 
初始 时 刻 ， 液 滴 随 机 分 布 在 x=3/8L 与 -=5/8L 之 
间 的 夹层 内 ， 之 后 液 滴 蒸 发 、 着 火 、 燃 烧 ， 并 沿 x 
轴 向 两 侧 扩散 。 本 文 所 选 时 刻 ， 液 滴 全 部 蒸发 ， 可 
认为 是 纯 气 相 燃 烧 。 详 细 信息 可 见 相 关 文献 上 9。 
本 文 使 用 的 化 学 反应 机 理 为 : 


110, + C,H, > 7CO, +8H,0 (8) 


Hh 


反应 速率 遵循 Arrhenius 定律 ， 其 参数 根据 
Westbrook" 的 研究 选择 ， 指 前 因子 取 


5.1x10" cm /mol.s， 活 化 能 取 30.0 Kcal/mol ， 庚 烷 


和 氧气 质量 分 数 的 指数 分 别 为 0.25 和 1.5。 具 体 的 
Arrhenius 表达 式 为 : 


Y 


Y m n 
w, = W v Aexp(-E/RT) p= p— 
| W W (9) 


fu 


m+n m n 
= Bp™KY,"Y, 


AP, A 为 指 前 因子 ，E 为 活化 能 ， 刺 为 组 分 
的 摩尔 质量 ，Y 为 组 分 的 化 学 计量 数 ，m 和 分 别 


为 组 分 的 指数 ，B = 4W "W v ，K =exp(-E/RT), 


图 1 为 所 选 时 刻 的 平均 温度 、 组 分 浓度 分 布 图 。 


图 1 平均 温度 、 组 分 浓度 分 布 曲线 


Figure 1. The distribution of the average temperature and 


mass fraction 

由 图 可 知 , 此 时 刻 , 中 间 区 域 为 高 温 的 反应 区 ， 
反应 物 C7Hi6 和 Os 均 大 量 存在 , 左右 两 侧 为 火焰 锋 
面 ， 是 新 鲜 气 体 与 已 燃气 体 的 分 界面 ， 火 焰 从 中 间 
向 两 侧 扩散 。 本 文 使 用 该 数据 库 ， 对 SOM 燃烧 模 
型 进行 先 验 性 检验 ， 所 采用 的 平均 数据 均 为 当地 统 
计 平 均值 , 即 yz 平 面 的 统计 平均 作为 x 处 的 当地 值 。 
3 SOM 燃烧 模型 检验 结 

本 文 根 据 二 阶 矩 思想 ,针对 表达 式 (2) 以 及 模 化 
方程 (4) 和 (7) 式 进行 了 检验 ， 以 下 是 检验 结果 。 
3. 1 反应 率 方程 的 检验 


表达 式 (2) 的 推导 采用 了 忽略 密度 脉动 和 三 阶 
脉动 关联 量 的 假设 ， 现 针对 该 假设 进行 检验 。 


0; 反 应 速率 /kg-m3.s1 


f=] 
[=] 
is} 


N 
[=] 
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xiL 
2 反应 率 方程 的 检验 
Figure 2. Validation ofthe mean reaction rate 
图 2 显示 了 分 别 忽略 密度 脉动 、 忽 略 脉 动 关 联 
量 以 及 二 者 同时 忽略 的 计算 值 与 与 准确 值 的 比较 结 
果 。 从 图 中 明显 的 看 出 ， 三 条 曲线 和 准确 值 吻合 的 


很 好 ， 其 中 最 大 的 误差 在 40% 之 内 。 由 此 可 知 ， 二 


行 的 ， 


阶 矩 方法 中 ， 忽 略 密度 脉动 和 三 阶 脉动 关联 量 是 可 
可 以 有 效 地 简化 模型 。 但 是 , 值得 注意 的 是 ， 


该 验证 所 采 | 


存在 密 


度 梯度 ， 当 使 用 党 


] 的 密度 是 当地 平均 密度 ， 在 x 方向 上 
密度 的 不 可 压 假设 时 ， 将 


联 量 输 运 方程 。 


会 对 结果 影响 很 大 ， 因 此 需要 


通过 对 精确 的 平均 反应 速率 的 展开 式 右 端 项 大 


E 导 可 压条 件 下 的 关 


关联 量 输 运 方程 ， 即 浓度 -浓度 关联 量 
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以 及 浓度 -反应 率 系 数 关联 量 方 程 (6) 中 的 精确 的 产 
生 项 、 扩 散 项 、 耗 散 项 以 及 化 学 反应 项 都 需要 进行 
检验 ， 图 3 显示 了 这 些 精 确 项 的 大 小 关系 。 


从 图 3(9) 可 以 看 出 ， 这 
量 级 上 ， 
产生 项 ， 最 后 是 耗 散 项 。 


小 的 比较 ， 发 现 平均 值 乘积 Ky"Y* 是 反应 速率 的 


主要 项 ， 决 定 了 反应 速率 的 方向 和 大 小 ; 


动 的 关联 量 ， 


rl 


火焰 面 处 则 存在 一 定 程度 的 峰 
对 于 三 阶 关 联 量 ， 除 了 在 火焰 


在 其 他 区 域 则 几乎 很 小 ， 


烧 模 型 中 ， 
的 ， 而 二 阶 脉动 关联 量 


fu ox 


而 含有 脉 


在 反应 区 有 相对 较 小 的 波动 ， 在 两 侧 


下 ， 是 不 可 忽略 的 ; 


用 处 有 少许 波动 外 ， 


可 以 忽略 。 
忽略 三 阶 脉动 关联 量 的 假设 是 可 以 接受 


因此 ，SOM 燃 
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3.2 模 化 方程 的 检验 


准确 的 模 化 。 


四 个 项 的 大 小 在 一 个 数 
其 中 ， 反 应 项 数值 最 大 ， 
图 3(b) 则 


其 次 是 扩散 项 、 
明显 地 显示 出 化 


学 反应 项 贡献 最 大 。 王 方 等 人 [的 不 可 压 DNS 的 
检验 结果 认为 ,化 学 反应 源 项 对 右 端 项 的 页 献 很 小 ， 


这 与 本 文 得 到 的 结果 有 所 差别 。 产 


生 该 结果 的 原因 


在 于 所 选取 的 数据 库 中 的 燃烧 热 释 率 不 同 ， 王 方 等 


人 的 反应 热 释 率 远 小 于 本 文 庚 烷 燃 烧 热 释 率 ， 二 者 


大 概 相 差 两 个 数量 级 ， 这 一 点 可 以 从 二 者 的 温度 曲 
即 王 方 等 人 "使 用 的 


线 看 出 ， 见 图 1 和 文献 [13]。 


槽 道 反应 流 温度 远 低 于 庚 烷 燃烧 温度 


王 方 等 人 统计 的 各 个 组 分 的 化 学 反应 速率 相对 较 


低 ， 低 佑 了 化 学 反应 源 项 的 页 献 。 


精确 项 /kgm3.s 


(b) 


3 ”关联 量 方程 右 端 项 大 小 检验 : (a) 浓 度 -浓度 脉动 关联 量 (b) 浓 度 -反应 率 系数 脉动 关联 量 (0O, 输 运 方程 ) 
Figure 3. Validation of the terms on the RHS of (a) Eq. (3) and (b) Eq. (6) for O2 


产生 项 /kg.m-3.s1 


扩散 项 /kg-m-3.s-1 


4 浓度 -浓度 模 化 项 检验 : (3) 产 生 项 (b) 扩 散 项 (0) 耗 散 项 (d) 反 应 项 


Figure 4. Validation of the terms (a) production term; (b) diffusion term; (c) dissipation term; (d) reaction term for Eq.(4) 


图 4、 图 5 分 别 显示 了 浓度 -浓度 脉动 关联 量 以 
及 浓度 -反应 率 系数 脉动 关联 量 中 模 化 项 的 检验 结 
Ro MB 4a, 4b, Sa 以 及 5b 可 以 看 出 ， 产 生 项 和 
扩散 项 的 模 化 效果 较 好 , 可 以 捕捉 到 精确 项 的 趋势 ， 


并且 可 以 在 两 侧 的 火焰 锋面 处 捕 提 到 两 个 峰 ， 其 中 
不 足 之 处 在 于 模 化 项 对 系数 的 选取 有 很 大 的 依赖 
性 ， 需 要 改进 。 图 4c 表示 了 浓度 -浓度 脉动 关联 量 
的 耗 散 项 的 模 化 结果 ， 从 图 中 可 以 明显 的 看 出 虽然 
模 化 结果 与 精确 项 的 基本 走势 相近 ， 但 是 精确 项 在 
火焰 锋面 处 存在 两 个 波峰 ， 而 模 化 项 只 有 一 个 峰 ， 
精确 项 两 个 波峰 的 产生 ， 是 由 于 火焰 面 两 侧 的 组 分 
梯度 脉动 引起 的 ， 而 使 用 组 分 脉动 和 注 动 能 的 耗 散 
来 封闭 , 是 难以 捕捉 这 个 双 峰 趋势 的 。 图 Se 表示 了 
浓度 -反应 率 系数 脉动 关联 量 的 耗 散 项 的 模 化 结 

与 图 4c 一 样 ， 结果 不 是 很 理想 , 但 是 由 于 其 在 输 运 


方程 中 的 贡献 较 小 ， 可 以 忽 


显示 了 浓度 -浓度 关联 量 


方程 的 反应 项 的 模 化 


可 以 明显 的 发 现 ， 化 学 反应 源 项 的 封闭 效果 不 是 很 
理想 。 原 因 在 于 反应 项 的 封闭 中 ， 没 有 考虑 反应 率 


各 其 带 来 的 误差 。 图 4d 


+ 四 
结果 ， 


系数 天 的 脉动 , 即 忽略 了 对 结果 影响 较 大 的 温度 脉 


动 。 因此 , 有 必要 对 化 学 反应 


硕 采取 新 的 封闭 方法 ， 


可 以 把 反应 率 系数 的 脉动 考虑 在 内 。 与 图 4d 相反 ， 


浓度 -反应 率 系数 关联 量 


好 ， 模 化 结果 可 以 在 中 
确 值 ， 但 在 两 侧 火 焰 锋 


闻 反 应 
处 却 


= 
= 


区 域 很 好 的 捕捉 
Minis, Ree 


的 反应 项 的 封闭 结果 比较 


到 精 
为 模 


拟 的 均匀 各 向 同性 清流 反应 中 ， 中 间 燃 烧 区 的 温度 


变化 不 会 很 大 ， 其 脉动 不 大 ， 
组 分 和 温度 的 急剧 变化 ， 存 在 较 大 的 梯度 ， 将 会 导 
致 反应 项 增 大 ， 因 此 需要 考虑 实际 组 分 、 温 度 梯 度 


Ew 
和 影响 。 


而 在 火焰 面 处 ， 


由 于 


5 浓度 -反应 率 系数 模 化 项 检验 (0; 输 运 方程 ): (a) 产 生 项 (D) 扩 散 项 〈c) 耗 散 项 (d) 反 应 项 


Figure 5. Validation of the terms (a) production term; (b) diffusion term; (c) dissipation term; (d) reaction term for Eq.(7) for O2 
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